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IMPROVED WAVEFRQNT CODING OPTICS 

CROSS-REFERENCE TO RELATED APPLICATIONS 

U.S. Patent No. 5,748,371, issued May 5, 1998 and entitled "Extended 
Depth of Field Optical Systems," and U.S. Patent No. 6,021,005, 
issued February 1, 2000 and entitled "Antialiasing Apparatus and 
Methods for Optical Imaging, "are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION: 

This invention relates to improved Wavefront Coding Optics for 
controlling focus related aberrations, and methods for designing 
such Wavefront Coding Optics. 

DESCRIPTION OF THE PRIOR ART : 

Wavefront Coding is a relatively new technique that is used to 
reduce the effects of misfocus in sampled imaging systems through 
the use of Wavefront Coding Optics which operate by applying 
aspheric phase variations to wavefronts of light from the object 
being imaged. Image processing of the resulting images is required 
in order to remove the spatial effects of the Wavefront Coding. The 
processed images are sharp and clear, as well as being relatively 
insensitive to the distance between the object and the detector. 
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Wavefront Coding is also used to control general focus related 
aberrations to enable simplified design of imaging systems as well 
as to provide anti-aliasing in sampled imaging systems. 

The Wavefront Coding Optics taught and described in prior art, 
beginning with U.S. Patent No. 5,748,371, issued May 5, 1998, were 
discovered by trial and error. The first operative Wavefront Coding 
mask applied a cubic phase function to the wavefront from the 
object. It was known that Wavefront Coding Optics, like the cubic 
mask, needed to apply aspheric, asymmetric phase variations to the 
wavefronts. 

Prior art related to Wavefront Coding systems includes a 
fundamental description of Wavefront Coding (U.S. Patent No. 
5,748,371), description of Wavefront Coding used for anti-aliasing 
(Anti-aliasing apparatus and methods for optical imaging, U.S. 
Patent #6,021,005, February 1, 2000), use of Wavefront Coding in 
projection systems (Apparatus and methods for extending depth of 
field in image projection systems, U.S. Patent #6,069,738, May 30, 
2000), and the combination of Wavefront Coding and amplitude 
apodizers (Apparatus and method for reducing imaging errors in 
imaging systems having an extended depth of field, U.S. Patent 
#6,097,856, August 1, 2000). 

The layout of a conventional Wavefront Coded imaging system is 
shown in Figure 1. Imaging Optics 104 collects light reflected or 
transmitted from Object 102. Wavefront Coding Optics 106 modify 
the phase of the light before detector 108. Wavefront Coding Optics 



106 comprise a cubic mask. Detector 108 can be analog film, CCD or 
CMOS detectors, etc. The image from detector 108 is spatially 
blurred because of Wavefront Coding Optics 106. Image processing 
110 is used to remove the spatial blur resulting in a final image, 
i.e., image processing 110 removes the Wavefront Coding applied by 
optics 106, thereby reversing the effects of optics 106, other than 
the increase in depth of field and depth of focus. The image before 
and after Image Processing 110 also is very insensitive to misfocus 
aberrations. These misfocus aberrations can be due to the Object 
102 being beyond the depth of field of the Imaging Optics 104, the 
detector 108 being beyond the depth of focus of the Imaging Optics 
104, or from Imaging Optics 104 having some combination of 
misfocus aberrations such as spherical aberration, chromatic 
aberration, petzval curvature, astigmatism, temperature or pressure 
related misfocus. 

Figure 2 describes a rectangularly separable prior art Wavefront 
Coding phase function that produces an extended depth of field. This 
phase function is a simple cubic phase function that is 
mathematically described, in normalized coordinates, as: 

cubic-phase(x.y) = 12 [ x 3 + y 3 ] 
lxl<1,lyl<1 

Other related forms of the cubic mask are described as: 



cubic-related-forms(x.y) 



= a[sign(x)lxl b + sign(y)lyl b ], 
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where 



Ixl < 1, lyl < 1, 



and 



sign(x) = +1 for x > 0, sign(x) = -1 otherwise 



These related forms trace out "cubic like" profiles of increasing 
slopes near the end of the aperture. 

The top plot of Figure 2 describes a 1D slice along an orthogonal axis 
of the cubic phase function. The lower plot of Figure 2 describes the 
contours of constant phase of this cubic phase function. 

Figure 3 shows MTFs as a function of misfocus for a system with no 
Wavefront Coding and for a system with the conventional Wavefront 
Coding cubic phase function of Figure 2. The normalized misfocus 
values are the same for both systems and are given as y = { 0 , 2 , 4}, 

where y = [ 2 pi W 20 ], and where W 20 is the conventional misfocus 

aberration coefficient in waves. MTFs with no Wavefront Coding 
(302) are seen to have a large change with misfocus. MTFs with the 
rectangularly separable cubic phase function (304) are seen to 
change much less with misfocus then the system with no Wavefront 



A non-separable prior art form of Wavefront Coding Optics, in 
normalized coordinates, is: 



Coding. 



non-separable-cubic-phase( p, q ) = p 3 cos( 3 q ) 



Ipl < 1, 0 < q < 2 pi 



This phase function has been shown to be useful for controlling 
misfocus and for minimizing optical power in high spatial 
frequencies, or antialiasing. When using a digital detector such as a 
CCD or CMOS device to capture image 108, optical power that is 
beyond the spatial frequency limit of the detector masquerades or 
"aliases" as low spatial frequency power. For example, say that the 
normalized spatial frequency limit of a digital detector is 0.5. As 
seen from Figure 3, the in-focus MTF from the conventional system 
with no Wavefront Coding can produce a considerable amount of 
optical power beyond this spatial frequency limit that can be 
aliased. By adding misfocus to the system without Wavefront Coding 
the amount of high spatial frequency optical power can be decreased, 
and aliasing reduced, as is well known. When using conventional 
Wavefront Coding, as shown in Figure 3, the amount of optical power 
that can be aliased can be decreased (304) compared to the system 
without Wavefront Coding (302). 

Image Processing function 110 essentially applies amplification and 
phase correction as a function of spatial frequency to restore the 
MTFs before processing to the in-focus MTF from the conventional 
system with no Wavefront Coding after processing, or to some other 
application specific MTF, if desired. In effect, the Image Processing 
function of Figure 1 removes the Wavefront Coding blur in the 
detected image. 

In practice the amplification applied by the Image Processing 



function increases the power of the deterministic image but also 
increases the power of the additive random noise as well. If Image 
Processing 110 is implemented as a linear digital filter then a 
useful measure of the increase of power of the additive random 
5 noise is called the Noise Gain of the digital filter. The concept of 

"noise gain" is commonly used in radar systems to describe the 
amount of noise power at the output of radar digital processors. 
Nonlinear implementations of Image Processing 110 have similar 
types of noise-related measures. The Noise Gain for a digital filter 

1 0 is defined as the ratio of the root-mean-square (RMS) value of the 
^ noise after filtering to the RMS value of the noise before filtering. 
J3 In general the Noise Gain is nearly always greater than one in 

£ Wavefront Coded systems. Assuming that the additive noise is 

yj uncorrelated white gaussian noise, the Noise Gain of a two 

1|§ dimensional linear digital filter can be shown to be equal to: 

■ H 

S Noise Gain = sqrt[ I Z f(i,k) 2 ] = sqrt[ I X IF(Wj,w k )l 2 ] 

□ where 

[ IX f(i.k) ] = F(0,0) = 1.0, 

f(i,k) is a spatial domain digital filter, F(Wj,w k ) is the equivalent 

2 0 frequency domain digital filter, and the first sum is over the index i 

or k and the second sum is over the other index. Indices (i,k) denote 
spatial domain coordinates while indices (Wj,w k ) denote frequency 
domain coordinates. The constraints that the sum of all values of 
the filter and the zero spatial frequency filter value both equal unity 



ensures that the zero spatial frequency components of the image 
(the background for example) are unchanged by the image processing. 

Wavefront Coded MTFs that have the highest values require the least 
amplification by the digital filter and hence the smallest Noise Gain. 

5 In practice the Wavefront Coding Optics that produce MTFs that have 

small changes over a desired amount of misfocus and also have the 
highest MTFs are considered the best and the most practical Optics 
for Wavefront Coding. Optics that produce MTFs that have small 
changes with misfocus but also very low MTFs are impractical due 
1 0 to very large Noise Gain of the resulting digital filters. Digital 

y filters with large Noise Gain will produce final images that have 

"fl unnecessarily high levels of noise. 

W While the conventional cubic Wavefront Coding mask does operate to 

fU increase depth of field and control focus related aberrations, there 

1 C§ remains a need in the art for improved Wavefront Codings Optics, 

hj which retain the capacity to reduce focus-related aberrations, while 

Q also producing high value MTFs. There also remains a need in the art 

i - 

^ for methods of designing such improved Wavefront Coding Optics. 



20 SUMMARY OF THE INVENTION 

An object of the present invention is to provide improved Wavefront 
Coding Optics, which retain the capacity to reduce focus-related 
aberrations, while producing high value MTFs, and to provide 
methods of designing such improved Wavefront Coding Optics. 
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Improved Wavefront Coding Optics according to the present 
invention, which apply a phase profile to the wavefront of light from 
an object to be imaged, retain their insensitivity to focus related 
aberrations, while increasing the resulting MTFs. Such improved 
Wavefront Coding Optics have the characteristic that the central 
portion of the applied phase profile is essentially constant, while 
the edges of the phase profile have alternating negative and positive 
phase regions (e.g. turn up and down at respective profile ends). 

To achieve higher MTFs, control misfocus and misfocus aberrations, 
and improve antialiasing characteristics, the central group of rays 
should be left unmodified by the Wavefront Coding Optics. In order 
to increase the light gathering (and possibly spatial resolution) of 
the full aperture system, the outer rays need to be modified. Only 
these outer rays need to be modified in order to increase the light 
gathering while keeping the depth of field and/or aliasing 
characteristics constant. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 (prior art) shows a conventional prior art Wavefront Coding 
imaging system. 

Figure 2 (prior art) shows a 1D plot of the Wavefront Coding cubic 
phase function of Figure 1 and a contour plot of the 2D 
representation of this function. 
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Figure 3 (prior art) shows MTFs as a function of misfocus for a 
system with no Wavefront Coding and with the conventional 
Wavefront Coding cubic phase function of Figure 2. 

Figure 4 (prior art) shows the ray paths for a system with no 
Wavefront Coding and with the conventional Wavefront Coding cubic 
phase function shown in Figure 2. 

Figure 5A shows the ray paths for a linear region phase function 
according to the present invention. 

Figure 5B shows the ray paths for a cubic region phase function 
according to the present invention. 

Figure 6 shows 1D plots of the linear region and cubic region phase 
profiles of Figures 5a and 5b, as well as contour plots of the 2D 
representations of these Wavefront Coding functions. 

Figure 7 shows MTFs as a function of misfocus for a system with no 
Wavefront Coding and with the linear region and cubic region Optics 
of Figures 5a and 5b. 

Figure 8 shows a contour plot of an improved non-separable phase 
function according to the present invention, a plot of a 1D slice 
through the non-separable phase function, and the MTFs as a function 
of misfocus for a non-separable cubic sector phase function. 

Figure 9 (prior art) shows a drawing of a conventional Cooke triplet 
lens. 
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Figures 10A and 10B (prior art) show contour plots of the exit pupils 
and corresponding MTFs as a function of field angle for the full 
aperture (10A) and stopped down (10B) Cooke triplet lens of Figure 
9. 

Figures 11A and 1 1 B show contour plots of the exit pupils and 
corresponding MTFs as a function of field angle for the prior art 
cubic phase Wavefront Coding system of Figure 1 (11 A) and the 
improved cubic region Wavefront Coding system of Figure 8 (1 1 B). 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

There are an infinite number of Wavefront Coding Optics that will 
reduce the variation in the resulting MTFs and PSFs of a given 
optical system due to misfocus or misfocus aberrations. Many of the 
possible optics are impractical in that the required Image 
Processing function 110 used to remove the Wavefront Coding blur 
from detected images would amplify the additive noise in practical 
images beyond an acceptable level. Improved forms of Wavefront 
Coding Optics that can control misfocus and misfocus aberrations, 
that can lead to higher MTFs, and that have improved antialiasing 
characteristics, as well as new methods of Wavefront Coding design, 
are shown in Figures 5 through 9. Use of these improved forms of 
Wavefront Coding Optics and design methods to control the misfocus 
aberrations with a Cooke triplet lens are shown in Figures 9 through 
11. 
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The improved Wavefront Coding Optics according to the present 
invention share the characteristic that the central region of the 
applied phase profile is essentially constant, while the edges of the 
phase profile have alternating positive and negative phase regions. 
5 Such Wavefront Coding Optics retain the ability to reduce focus 

related aberrations, but also have considerably higher MTFs than 
traditional Wavefront Coding Optics, thus reducing noise in the final 
images produced. 

Wavefront Coding Optics (in the form of aspheric optical elements) 

1 0 are placed at or near the aperture stop of optical systems (or at or 
□ near an image of the aperture stop) in order to redirect light rays as 
S a function of spatial position. The aspheric optical elements can be 
ry constructed of optical glass or plastic with varying thickness 

jj and/or index of refraction. The optics can also be implemented with 

is 3 

t3 shaped mirrors, spatial light modulators, holograms, or micro mirror 

S devices. U.S. Patent No. 6,021,005, issued February 1, 2000 and 

g entitled "Antialiasing Apparatus and Methods for Optical Imaging," 

p provides descriptions of a variety of devices for applying variations 

to a wavefront of light from an object. 

20 Light rays from an ideal thin lens with no Wavefront Coding 

converging to focus 50mm from the lens are shown on the top graph 
of Figure 4. All rays from the ideal system without Wavefront 
Coding travel towards the best focus position on the optical axis. 
The light rays from a conventional (prior art) rectangularly 

2 5 separable cubic phase system are shown on the bottom graph of 
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Figure 4. Notice that the rays from the top half of this lens cross 
the optical axis beyond the best focus point of the conventional lens 
(or 50mm). The rays from the bottom half of this lens cross the 
optical axis before the best focus point of the conventional lens. 

Since no two rays from the prior art Wavefront Coding system cross 
the optical axis at the same point, every ray of the prior art 
Wavefront Coding cubic phase system is being modified when 
compared to the system not using Wavefront Coding (except the on- 
axis zero-slope ray). To achieve higher MTFs, control misfocus and 
misfocus aberrations, and improve antialiasing characteristics, the 
central group of rays should be left unmodified. 

Consider a full aperture and a stopped down conventional optical 
system that does not use Wavefront Coding. Assume that a 
particular application has a depth of field (or depth of focus) and/or 
antialiasing requirement that the full aperture system cannot meet. 
It is well known that stopping down the lens will increase the depth 
of field of the system while reducing the potential spatial 
resolution of the optics. Stopping down the aperture also reduces 
the optical power available at the detector. There is a particular 
stopped down aperture where the imaging system has a best fit to 
the depth of field and/or antialiasing characteristics required. The 
light rays that pass through the aperture of the stopped down 
system are then considered suitable from a depth of field and/or 
antialiasing perspective for the particular application. 

In order to increase the light gathering (and possibly spatial 
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resolution) of the full aperture system, the rays of the full aperture 
system that lay outside the stopped down aperture need to be 
modified. Only those rays that lay outside the stopped down 
aperture need to be modified in order to increase the light gathering 
while keeping the depth of field and/or aliasing characteristics 
constant. When the rays within the stopped down aperture are 
modified, as in all prior art Wavefront Coding Optics, the resulting 
MTFs are not as high as possible, the resulting Noise Gains are not as 
low as possible, and the resulting images are more noisy than 
necessary. 

By not modifying the central rays of Wavefront Coding systems the 
resulting MTFs below the spatial frequency limit of the digital 
detector can be increased over prior art Wavefront Coding systems. 
Only the rays outside of the central region of the aperture need to be 
modified in order to control misfocus or misfocus aberrations. The 
central region is defined as the general region of a stopped down 
aperture where a system would have a suitable depth of field, depth 
of focus, or antialiasing characteristics for a particular application, 
albeit with reduced light gathering and spatial resolution. 

Improved rectangularly separable Wavefront Coding Optics can be 
most generally described mathematically through phase functions, in 
normalized coordinates, that have the form: 

phase(x,y) = Z[ U(lxl/A xj ) G xj (x) + U(lyl/A yi ) G yj (y) ] 

where 
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Ixl < 1, lyl < 1 
i = 1, 2, ... , N 

and where 

U(z) = 1 if z >1 , U(z) = 0 otherwise 
5 0< A xj <1, 0< A yj <1 

The sum is over the index i. The function U(lxl/A x ) is a zero/one 
step function that has a value of zero inside of an aperture of length 
2 A x and has a value of one outside of this aperture. The functions 
3 G x and G y are general functions that modify the rays of the system 

;,: 5 
tssr. 

1J9 outside of the specific aperture defined by A x and A The aperture 

yj shape in this form is described by a rectangle for mathematical 

ry convenience, but in general can be described by any closed shape. 

a 

□ For example, instead of a rectangular aperture, a circular, elliptical, 

y or multi-sided polygonal aperture can also be be used. 

S : 
fl 

13 In order for the phase function to control misfocus effects, the 

r"*~ 

phase functions G x and G y should be designed so that groups of rays 
from specific regions of the aperture cross the optical axis either 
before or after the point of best focus when no Wavefront Coding is 
used. 

20 With these concepts, a number of improved rectangularly separable 

Wavefront Coding Optics can be formed. These optics differ in the 
composition of the general functions G x and G y . For example, a 
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linear phase region optical system can be described as: 

linear-phase-region(x,y) = U(lxl/A x ) sign(x)(lxl-A x )/(1 -A x ) 



+ U(lyl/A y ) sign(y)(lyl-A y )/(1-A y ) 

where 

Ixl < 1, lyl < 1 
U(z) = 1 if z <1, U(z) = 0 otherwise 
sign(z) = +1 for z > 0, sign(z) = -1 otherwise 
0< A x <1, 0< A y <1 

The linear phase region system has zero phase inside of the aperture 
defined by A x and A y with linearly changing phase as a function of 

spatial position variables x and y outside of this aperture. The 
linear phase region system delivers extended depth of field with 
high MTFs below the spatial frequency limit of the digital detector. 
This system also has a very simple physical form that can be 
preferred over smoothly varying forms because of the fabrication 
process or physical implementation being used. 

For example, if fabrication of the Wavefront Coding Optics is done 
with a precision milling machine, the linear phase region optics 
would be preferred over the cubic phase optics because the linear 
phase region optics have only fixed surface slopes while the cubic 
phase optics have continuously changing surface slopes. Often, 

1 5 



optical elements with fixed surface slopes are easier to fabricate 
than elements with continuously changing slopes. If the Wavefront 
Coding optical surface is implemented with micro mirrors or 
similar, the small number of surface slopes of the linear phase 
region system can be more easily implemented than the continuously 
changing surface slopes of the prior art cubic phase optics. 

Another version of improved rectangularly separable Wavefront 
Coding Optics is mathematically described by the sum of powers 
region phase function. This phase function is defined as: 

sum-of-powers-region(x,y) = 

Z U(lxl/A xi ) a, sign(x) [ (Ixl -A xi )/(1 -A xj ) ] Bi 

+ 2 U(lyl/A yj ) Xi sign(y) [ (lyl -A yi )/(1 -A yi ) ] 5i 

where 

Ixl < 1, lyl < 1 
1 = 1,2, ... , N 
U(z) = 1 if z < 1, U(z) = 0 otherwise 
sign(z) = +1 for z > 0, sign(z) = -1 otherwise 
0< A xj <1, 0< A yj <1 

and where the summations are over the index i. 
Examples of the linear phase region systems and the sum of powers 
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region systems are found in Figures 5, 6, and 7. The graph of Figure 
5A shows the ray paths from one dimension of the linear phase 
region system. The graph of Figure 5B shows the ray paths for one 
dimension of a cubic phase region system. A cubic phase region 
system is sum of powers phase region system with only a single 
term and the exponent on this term having the value of 3, or 3 } = 8j = 

3, and otj = % { = 0 for i * 1 . 

From Figure 5A we see that the rays from the central region of the 
lens (with this central region extending from -2.5 to +2. 5mm) are 
not modified by the linear phase region system. All the unmodified 
rays travel towards the optical axis at the best focus position 
50mm to the right of the lens. Only the outer rays from the lens, 
those beyond +/-2.5mm, are modified by the linear phase region 
function. The rays from the top of the lens are modified so that they 
cross the optical axis at points beyond the best focus position. The 
rays from the bottom of the lens cross the optical axis at points 
before the best focus point. 

From Figure 5B, the cubic phase region system also does not modify 
the central region (extending from -2.5mm to +2. 5mm) of the lens. 
The remaining rays are modified in a spatially varying manner such 
that one set of rays (from the top of the lens) crosses the optical 
axis at points beyond the best focus point and another set of rays 
(from the bottom of the lens) crosses the optical axis at points 
before the best focus position. 
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The linear phase region system and the sum of powers phase region 
system can consist of different regions that have different slopes 
and/or different number of terms. For example, the linear phase 
region system, as shown in Figure 5A, instead of having a zero phase 
central region and two linearly changing regions, could have a zero 
phase central region and more than two linearly changing regions. 
Additional linear regions could direct the rays to different positions 
in order to better control misfocus, misfocus aberrations, and 
antialiasing characteristics. The size and shape of the misfocus 
PSFs can also be visualized and controlled by controlling regions of 
phase function and the corresponding rays. When only two 
asymmetric regions of the phase function are being used (as in 
Figure 5A), manual optimization is possible. When the number of 
regions exceeds two, then the number, size, and phase of each region 
is generally best determined by computer optimization. 

Figure 6 gives another view of the linear phase region and cubic 
phase region Wavefront Coding systems. The top graph of Figure 6 
shows a 1D slice of phase functions describing the optics of both 
systems along one of the orthogonal axes. The 1D form of this linear 
phase region system is: 

3.86 U(lxl/A x ) sign(x) [ (Ixl - A x )/(1- A x ) ] 

where 

A x = 2/3 
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Ixl < 1 



U(z) = 1 if z < 1, U(z) = 0 otherwise 
sign(z) = +1 for z > 0, sign(z) = -1 otherwise 
The 1D form for the cubic phase region system is specifically: 
5 U(lxl/A x ) sign(x) [ (Ixl -A x )/(1-A x ) ] 3 
A x = 1/3 
Ixl < 1 

U(z) = 1 if z < 1, U(z) = 0 otherwise 

sign(z) = +1 for z > 0, sign(z) = -1 otherwise 

The 2D contours of constant phase for these phase functions are 
shown in the bottom of Figure 6. These contours clearly show that 
the phase within the central regions of these systems are constant; 
i.e. the central rays of the corresponding Wavefront Coding systems 
are not modified. The phase near the edges of the linear phase 
region system increases/decreases linearly while the phase near the 
edges of the cubic phase region system increases/decreases as a 
cubic function. 

Figure 7 shows the misfocus MTFs as a function of normalized 
misfocus for the system with no Wavefront Coding, the linear phase 
region system, and the cubic phase region system. The normalized 
misfocus values are y = { 0 , 2 , 4 }, as used in Figure 3. The MTFs 
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from both the linear phase region and cubic phase region systems are 
very close to each other (and thus are not specifically distinguished 
in Figure 7) and are very insensitive to misfocus, especially when 
compared to the system using no Wavefront Coding. Compare the 
heights of the improved Wavefront Coded MTFs of Figure 7 to those 
of the prior art cubic phase Wavefront Coded MTFs shown in Figure 3. 
The MTFs from both the linear phase region and the cubic phase 
region systems have higher MTFs than the conventional cubic phase 
function MTFs while all the Wavefront Coding MTFs are essentially 
insensitive to misfocus. In terms of antialiasing, the MTFs from the 
improved linear phase region and cubic phase region systems have 
higher values of MTFs, compared to the prior art cubic phase MTFs, 
for the lower spatial frequencies that typically are not aliased, 
while also reducing the high spatial frequency optical power when 
compared to the system with no Wavefront Coding. 

Improved non-separable Wavefront Coded Optics can also be 
described. General non-separable Wavefront Coded Optics can be 
mathematically defined through the form: 



phase(p,0) = E CKp/Qj) Gj( p,9 ) 



where 



Ipl < 1, 0 < 0 < 2 pi 



1 = 1,2, ... N 
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where Q(z) = 1 if z > 1 , Q(z) = 0 otherwise 



0 < Qj < 1 



and where the sum is over the index i. The function Q(p/Q) is a 
zero/one function that allows the central rays from a region with 
radius less than Q to be unchanged. The function Gj( p,0) is a general 
phase function that has regions of positive and negative phase so 
that regions of rays are made to either cross the optical axis before 
the best focus image point, or after the best focus image point. 

One general example of improved non-separable Wavefront Coded 
Optics is mathematically described as : 

non-separable-phase(p,G) = Z Q(p/^) otj p m cos( Wj 9 - ) 

where 



Another example of improved non-separable Wavefront Coded Optics 
is the non-separable sum of powers form given by: 



Ipl < 1, 0 < 9 < 2 pi 



i = 1,2, ... N 



and where Q(z) = 1 if z > 1, Q(z) = 0 otherwise 



0 < Qj < 1 



non-separable-sum-of-powers(p,6) = 



21 



Z Q(p/QJaj sign 



angle 



(M9 - offset) [(p -a^)/(i-aj] m 



where 



Ipl < 1, 0 < 0 < 2 pi 



1 = 1,2, ... ,N 



and where Q(z) = 1 if z > 1, Q(z) = 0 otherwise 



0 < < 1 



sign angle ( O ) = +1 if 0 < O < pi, sign angle ( O ) = -1 otherwise 

where the integer M controls the number of +/- sectors used and 
where the sum is over the variable L 

An example of the improved non-separable sum of powers Wavefront 
Coded Optics is given in Figure 8. The upper left graph in Figure 8 is 
a contour plot of constant phase of the non-separable cubic phase 
sector system. The cubic phase sector system is the non-separable 
sum-of-powers system with one exponential term, with 3 = 3. M = 
1, and offset = pi/2. The parameters of this non-separable cubic 
phase sector system are: 



4.48 Q(p/Q) sign 



angle 



(e - pi/2) [( P -a)/(i-a)]3 



a = 0.2 



Ipl < 1, 0 < e < 2 pi 



and where 

Q(z) = 1 if z > 1 , Q(z) = 0 otherwise 

s '9 n angle( * ) = +1 if 0 < * < P'. si 9 n angle( ° ) = " 1 Otherwise 

The top right graph of Figure 8 shows a 1D slice through the cubic 
phase sector system. The bottom graph of Figure 8 shows the 
misfocus MTFs of a system with no Wavefront Coding and with the 
non-separable cubic phase sector system. Again the normalized 
misfocus values are \|/ = { 0 , 2 , 4 }, as used in Figures 3 and 7. 
These MTFs are 1D slices from MTFs formed with circular apertures. 
The MTFs for the system with the non-separable cubic phase sector 
system are seen to be very insensitive to misfocus effects, 
especially when compared to the MTFs from the system with no 
Wavefront Coding. 

Figures 9, 10, and 11 describe one example of using the improved 
Wavefront Coding Optics to control misfocus aberrations. This 
example shows control of field-dependent aberrations with a 
common Cooke triplet lens used for visible light. See Modern Optical 
Engineering (Warren J. Smith, McGraw-Hill, Inc, NY, 1990) for more 
information on the Cooke triplet. A drawing of this lens is given in 
Figure 9. The Noise Gain resulting from the improved Wavefront 
Coding Optics, as compared to prior art Wavefront Coding Optics, are 
shown to be a factor of 2.5 less than prior art systems with this 
lens. This reduced noise gain directly translates into final images 
that have 2.5 times less noise that prior art systems. 



This triplet uses all spherical surfaces except for the second 
surface of the second or middle lens, which contains the Wavefront 
Coding Optics. The prescription of the triplet lens without 
Wavefront Coding Optics is given by: 

Surface Radius Thickness Glass Diameter 



# 



Object 




Infinity 






1 


22.01359 


3.258956 


SK16 


17.10986 


2 


-435.7604 


6.007551 




15.98699 


3 


-22.21328 


0.9999746 


F2 


7.990108 


4 * 


20.29192 


4.750409 




7.731661 


5 


79.6836 


2.952076 


SK16 


13.78242 


6 


-18.3783 


42.18171 




14.5162 


Image 











* Surface #4 is the aperture stop 

All dimensions are given in units of mm. The focal length of this 
lens is 50mm, the full aperture F/# is 5, the half field of view is 20 
degrees. This lens is used in a system with a digital gray scale 
detector. The pixel pitch is 7.6 microns with 100% fill factor 
square pixels. The spatial frequency limit of this detector is 65.8 
Ip/mm. The ambient temperature is considered fixed at 20 degrees 
C with one atmosphere of pressure. Ideal fabrication and assembly 
are also assumed. 



Even though the lens of Figure 9 is a fairly complicated multi-glass 
optical system, this system suffers from the monochromatic 
misfocus aberrations of spherical aberration, petzval curvature, and 
astigmatism that limit off-axis performance. If the lens were 
designed with a single optical material, instead of the two different 
materials actually used, the lens would also suffer increased 
chromatic aberration. If plastic optics were used in place of glass 
optics the lens would also suffer from temperature related misfocus 
effects. 

Figure 10A (prior art) shows the performance of the system of 
Figure 9 operating at F/5 (wide open) and Figure 10B shows the 
system of Figure 9 with the aperture stopped down to F/19. Each of 
the following graphs are related to measures of the lens with green 
light. A contour plot of the on-axis exit pupil optical path 
difference (OPD) for the F/5 system is given on the left of Figure 
10A. The peak-to-valley OPD for this exit pupil has a value of 
approximately 0.7 X. The right plot of Figure 10A describes the MTFs 

for the F/5 system as a function of field angle for on-axis, 14 
degrees, and 20 degrees. Included in these MTFs, and all following 
MTFs, is the pixel MTF related to the 100% fill factor 7.6 micron 
square pixel. Only 1D slices along the horizontal axis of the 2D 
diffraction limited MTFs and the 2D Wavefront Coded MTFs are 
shown in Figures 10 and 11. The MTFs as a function of field angle 
are seen to vary greatly due to the presence of aberrations in the 
conventional F/5 system. 



The effect of stopping down the lens from F/5 to F/19 is seen in the 
graphs of Figure 10B. The peak-to-valley OPD of the on-axis stopped 
down exit pupil is reduced to 0.12 X. The MTFs as a function of field 
angle of the stopped down system are seen to change very little, 
5 compared to the full aperture system. Stopping down the aperture is 

one way of controlling the aberrations of this system, although at a 
reduction factor of (5/1 9) 2 or a loss of 93% of the optical power 
that was captured by the full aperture system. 

Figures 11 A and 1 1 B show the performance of the Cooke triplet of 
1 0 Figure 9, modified to utilize Wavefront Coding. The Wavefront 

y Coding Optics for this example were modelled as being added to the 

£J aperture stop of the system which is at second surface of the second 

[V element, or surface #4. The surface height equations for the prior 

J3 art Wavefront Coding cubic phase system used in Figure 11 A are: 

1g Z prior art (x,y) = 8.6061 E-5 {Ixl^ + lyl 3 } 

Ixl < 5.45, lyl < 5.45 



where the surface height Z prjor art (x,y) is given in mm. The optical 
area that is used is a circle of radius 3.86 mm. A square aperture 
can also be used in practice. 

20 The surface equations for the improved Wavefront Coding cubic 

phase region system used in Figure 1 1 B are: 



Z improved< x 'y) = 7.6535E-5 U(lxl) sign(x) [1x1-1]-= 
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+ 7.653 E-5 U(lyl) sign(y) [lyl-1] 3 



Ixl < 5.45, lyl < 5.45 
U(z) = 1 if z > 1, U(z) = 0 otherwise 
sign(z) = +1 for z > 0, sign(z) = -1 otherwise 

where again the surface height Z improved (x,y) is given in mm, the 

optical area that is used is a circle of radius 3.86 mm, and a square 
aperture can also be used in practice. 

The graphs of Figure 11A describe use of the prior art rectangularly 
separable Wavefront Coding cubic phase Optics. The graphs of Figure 
1 1 B describe use of the improved rectangularly separable cubic 
phase region Wavefront Coding Optics. Use of the improved 
Wavefront Coding Optics increases the height of the MTFs before 
image processing 110 and therefore drastically reduces the Noise 
Gain of the digital filters needed to give the ideal performance after 
processing, compared to the prior art optics. The improved MTFs 
also show increased antialiasing performance by having higher 
optical power below the detector spatial frequency cutoff, with 
greatly reduced MTFs beyond the detector cutoff, when compared to 
the conventional full aperture system. 

The prior art cubic phase system was designed to have a peak-to- 
valley OPD over the on-axis wide open (F/5) exit pupil of 
approximately 9.5 X. The resulting MTFs, before image processing, 



as a function of field angle are essentially constant, as shown in 
Figure 11 A, especially compared to the full aperture F/5 system 
without Wavefront Coding of Figure 10A. The diffraction limited 
MTF for this system was chosen as the desired system performance 
of the Wavefront Coding systems after Image Processing 110. In 
general, after image processing the resulting PSFs and MTFs in the 
Wavefront Coding system can have nearly any form. The Image 
Processing function 110, for this example, implements a 2D linear 
filter to transform the Wavefront Coding MTFs before image 
processing to MTFs after image processing that closely match the 
ideal diffraction limited MTF below the detector spatial frequency 
cutoff of 65 Ip/mm. The resulting Noise Gain values of the digital 
filters are then used as a figure of merit to judge the two 
dimensional height of the MTFs compared to the desired diffraction- 
limited MTFs. For the prior art cubic phase system of Figure 11A the 
Noise Gain of the resulting 2D digital filter is 8.1. 

The graphs of Figure 11 B describe use of the improved cubic phase 
region Wavefront Coding Optics. Since the stopped down F/19 
system of Figure 10B has suitable performance for rays within a 2 
mm diameter aperture, the cubic phase region system is constant, or 
has zero phase, over the square aperture region of +/- 1.0 mm. This 
zero phase region corresponds to the aperture of the stopped down 
F/19 system. 

The zero phase region could just as easily be formed into a circle or 
other geometric shape depending on the application and processing 
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to be used. The square aperture is more consistent with the 
rectangularly separable nature of the cubic phase region system than 
a non-separable circular region. 

The parameters of the cubic phase region system of Figure 1 1 B were 
designed to have similar insensitivity to off-axis misfocus 
aberrations as the prior art cubic phase system of Figure 11 A. This 
results in a peak-to-valley on-axis OPD also approximately 9.5 X. 

The contour plot of the exit pupil on the left of Figure 1 1 B clearly 
shows a large region near the center of the aperture that has zero 
phase as compared to the contour plot from the prior art cubic phase 
system that is not optically constant over any region in the aperture. 
The MTFs, before image processing, as a function of field angle for 
the improved cubic phase region system are seen to be essentially 
constant. The height of the MTFs resulting from the improved cubic 
phase region system is also much higher than those from the prior 
art cubic phase system below the spatial frequency band limit of the 
7.6 micron detector, or 65 Ip/mm. 

The Noise Gain of the 2D digital filter needed to match the 
performance of the MTFs after filtering to that of the diffraction- 
limited system has a value of approximately 3.2 for the improved 
cubic phase region system. Thus, the improved cubic phase region 
system produces nearly ideal performance in the control of the field 
dependent aberrations and also drastically reduces the digital filter 
Noise Gain from the prior art system by a factor of (8.1/3.2), or 
approximately 2.5. Thus, the additive noise power in the final 



images after Image Processing 110 will be 2.5 times larger with the 
prior art system than with the improved cubic phase region 
Wavefront Coding system. 

This large decrease in Noise Gain will result in substantially less 
noisy final images from the Improved Wavefront Coding Optics as 
compared to prior art optics. Or, for a constant amount of Noise 
Gain, the improved Wavefront Coding Optics can control a much 
larger degree of misfocus than can the prior art optics. 

Although not shown, th|s lens system when modified with Wavefront 
Coding also substarf(ictlly corrects misfocus arising from chromatic 
effects and temperature related effects, as well as reducing the 
tolerance of the system to manufacture and assembly errors. 

What is claimed is: 
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